CERN-TH.7540/94 
SHEP 95-10 



hep-ph/95042r7 



The Invisible Renormalon 



Guido Martinelli'^'*, Matthias Neuberf^ and Chris T. Sachrajda 

" Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

^ Department of Physics, University of Southampton 
Southampton S017 IB J, United Kingdom 

Abstract 

We study the structure of renormalons in the Heavy Quark Effective 
Theory, by expanding the heavy quark propagator in powers of 1/ uiq . We 
demonstrate that the way in which renormalons appear depends on the 
regularisation scheme used to define the effective theory. In order to inves- 
tigate the relation between ultraviolet renormalons and power divergences 
of matrix elements of higher-dimensional operators in the heavy quark 
expansion, we perform calculations in dimensional regularisation and in 
three different cut-off regularisation schemes. In the case of the kinetic 
energy operator, we find that the leading ultraviolet renormalon which 
corresponds to a quadratic divergence, is absent in all but one (the lattice) 
regularisation scheme. The nature of this "invisible renormalon" remains 
unclear. 
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1 Introduction 



Several recent papers have been devoted to the study of renormalons in the Heavy 
Quark Effective Theory (HQET) |T|-||§]. These studies exploit and extend the 
understanding of the relevance of renormalons in field theory, and for Operator 
Product Expansions (OPE) in particular, developed in refs. The present 



interest was initiated by the observation that the pole mass of a heavy quark, 
which is an important parameter in the heavy quark expansion, has an intrinsic 
ambiguity of order Aqcd due to the presence of infrared (IR) renormalons [|I], 
1^. In addition to the infrared renormalons and other non-perturbative effects 
present in QCD, renormalon ambiguities also arise in the HQET as a result of the 
expansion in powers of I/uiq. A consequence of these additional renormalons is 
the appearance of non-perturbative effects in the Wilson coefficients which relate 
operators in the effective theory to operators in the full theory (QCD). This 
problem is common to all effective theories obtained from an expansion in inverse 
powers of a large scale. Since these effects are an artefact of the construction of 
the effective theory and are absent in the original theory, they have to cancel 
in predictions for physical quantities. In the context of the HQET, it has been 
established that the IR renormalon ambiguities in the Wilson coefficients are 
cancelled by ultraviolet (UV) renormalon ambiguities in the matrix elements 
of the operators in the effective theory. How this cancellation occurs was traced 
explicitly, to order l/mg in the heavy quark expansion, for both inclusive f^, |^, |1 
and exclusive weak decays of hadrons containing a heavy quark. A non- 
perturbative approach for the elimination of the ambiguities in the definition of 
the pole mass, or more generally in the matrix elements of the operators appearing 



in the HQET, has been proposed in ref. |15 



In the present paper we study in more detail the structure of renormalons in 
the HQET. To be specific, we consider QCD Green functions that depend on a 
large scale tuq (the heavy quark mass) and a small scale k (the so-called resid- 
ual momentum). The purpose of the heavy quark expansion is to disentangle 
the physics on different length scales by introducing a factorisation scale A such 
that k <^ X rriQ. Contributions from virtual momenta above A are calculable 
in perturbation theory and attributed to Wilson coefficients, whereas contribu- 
tions from virtual momenta below A are contained in the matrix elements of the 
operators in the effective theory. If this program is performed with a "hard" 
factorisation scale, these matrix elements will, for dimensional reasons, diverge 
as powers of the UV cut-off A. The appearance of power divergences leads to 
non-perturbative ambiguities, since factors such as 



give non- vanishing contributions, which do not appear in perturbation theory [|6 



For practical reasons, however, one usually calculates the Wilson coefficients us- 
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ing dimensional regularisation. In this case, by definition, power divergences do 
not appear, and the "hard" factorisation scale A is replaced by a "soft" renormal- 
isation scale /i. In such a scheme it is unavoidable that the Wilson coefficients 
receive contributions from momenta below /i (so-called IR renormalons) , and 
the matrix elements receive contributions from momenta above fi (so-called UV 
renormalons) . These contributions lead to a factorial growth of the coefficients in 
the perturbative expansion of the Wilson coefficients and matrix elements. The 
corresponding perturbative series are divergent and not Borel summable. Since 
the choice of the renormalisation scale is arbitrary, however, the effects of IR and 
UV renormalons must cancel each other if one combines the perturbative series 
for the coefficient functions and matrix elements. 

Our main point is to demonstrate that the way in which renormalons appear 
in the HQET (and, by the same argument, in other effective field theories) is not 
universal, but depends on the regularisation scheme used to define the effective 
theory. To this end, we investigate in detail the relation between UV renormalons 
and power divergences of hadronic matrix elements of higher-dimensional opera- 
tors. We study explicitly the l/mg expansion of the heavy quark propagator in 
different regularisation schemes: dimensional regularisation and three different 
schemes with a hard UV cut-off (Pauli-Villars, momentum flow [0] and lattice 
regularisation). 

In general, from the degree of divergence of the matrix elements of an opera- 
tor in the HQET one can deduce the position of the UV renormalon singularities 
in the Borel transform of these matrix elements (see section 2 for a detailed 
explanation). By dimensional arguments we expect that the kinetic energy op- 
erator for a heavy quark is quadratically divergent, since it can mix with the 
identity operator under renormalisation. Indeed, this mixing has been observed 
in lattice perturbation theory ||T6[|. The corresponding renormalon singularity is 
absent, however.]^ This "missing" singularity is what we have called the "in- 
visible renormalon". Even more surprising is the absence, at one-loop order, of 
the quadratic divergence in the two other schemes considered in our study, i.e. 
the Pauli-Villars and momentum flow regularisations. We have been unable to 
understand whether there is a symmetry, broken by the lattice regularisation, 
which prevents the appearance of the quadratic divergence in the other cases, or 
whether this is an artefact of one- loop perturbation theory. In the latter case, 
the quadratic divergence would appear at higher orders. Given the relevance of 
the kinetic energy operator for heavy hadron spectroscopy and inclusive decays 



^-[^, an understanding of this puzzle is very important. 
The remainder of this paper is organised as follows: In sect. we summarise 
and discuss the main results of our study. Technical details and explicit calcula- 
tions in different regularisation schemes are presented in sects. |^-^. In sect, ^we 

^The absence of the corresponding IR renormalon in the pole mass was previously noted in 
ref. §. 



2 



give our conclusions. 



2 Renormalons and power divergences 

Perturbative expansions in QCD are asymptotic and, in general, not Borel sum- 
mable. When one tries to resum a perturbative series to all orders, one encounters 
ambiguities, which indicate that perturbation theory is by itself incomplete and 
must be supplemented by non-perturbative corrections. A convenient way to 
analyse these ambiguities is to consider the Borel transform S{u) of a series 
S{as{ft)) with respect to the coupling constant 0. Formally, the Borel sum of 
the series can be defined by the integral^ 

^(a.(/.))= /d«exp(-^^)^H, (2) 







where Po = 11 — -^Uf is the first coefficient of the /^-function. However, if the 
Borel transform contains singularities on the integration contour, the result of 
the integration depends on the regularisation prescription, and S{as{fj,)) is not 
uniquely defined in terms of S{u). In QCD, one source of such singularities 
are the higher-order diagrams in which a virtual gluon line with momentum k 
is dressed by a number of fermion, gluon and ghost loops. More precisely, one 
has to consider a gauge-invariant generalisation of such diagrams, which is called 
a renormalon chain. Effectively, this introduces the running coupling constant 
as{k) at the vertices. When as{k) is expressed in terms of the coupling constant 
renormalised at a fixed scale /i, the appearance of powers of large logarithms 
leads to a factorial divergence in the expansion coefficients of the perturbative 
series. Associated with this are renormalon singularities in the Borel transform 
S{u). Depending on whether these singularities are related to the region of large 
or small virtual momenta, they are referred to as UV or IR renormalons p|]-[P^. 



Although the resummation of renormalon chains corresponds to only a par- 
tial resummation of the perturbative series, it elucidates many non-perturbative 
effects and thus provides an interesting, non-trivial approximation. If one ac- 
cepts this so-called "bubble approximation", the renormalon singularities oc- 
cur as poles on the real axis in the Borel plane.0 Let us denote the position 
of the nearest pole on the positive w-axis by mq, i.e. S{u) = ro/{u ~ uq) + 
terms that are regular at u = uq. Then a measure of the ambiguity in the Borel 
integral (0) is given by 

AC f 47rao \ /Aqcd\'"° 



^We follow the notations and definitions of ref. 

•^When the calculations are extended beyond the "bubble approximation" , the poles become 
replaced by branch points of cut singularities. 
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In the last step we have used the one-loop expression for the running couphng 
constant. 

As mentioned above, in dimensional regularisation power divergences are 
hidden because of the absence of an intrinsic mass scale in the computation. 
However, perturbation theory "knows" about these divergences in the form of 
renormalon singularities in the Borel plane. In fact, there is a one-to-one cor- 
respondence between the structure of the UV renormalon poles in dimensional 
regularisation and the power divergences in regularisation schemes with a hard 
cut-off. For instance, if a quantity is linearly divergent in one-loop perturbation 
theory, then its Borel transform in the bubble approximation is logarithmically 
divergent at m = 1/2, corresponding to a pole singularity at this point. A similar 
argument is expected to hold for quadratic or higher-order power divergences. To 
illustrate this point we present results for the renormalised inverse heavy quark 
propagator S~^{k) in the HQET. Here k = pq —rnqv is the residual momentum, 
which stays finite in the limit niq oo, and f is a four- velocity vector (f^ = 1), 
which is usually identified with the velocity of the hadron containing the heavy 
quark. In dimensional regularisation, renormalon effects appear as singularities 
in the Borel transform S~^{u). In the bubble approximation one findsQ 



~ fir* 
S-i{k,u) = vk\6{u)^ 



+ 7; Rziu] 



/i / r(2 + M) 2u 



(4) 

in the Landau gauge. The scheme- dependent function Rz{u), which for renormal- 
isation schemes with analytic counterterms (such as MS-like schemes) is entire 
in the complex u-plane, is irrelevant for our purposes. Expanding the above ex- 
pression around u = and inserting the result into the Borel integral (|]), one 
finds 

S-^\k) = ^; . |l + ^ In + const.) + 0(a^)| , (5) 

where Cp = 4/3 is a colour factor. The F-functions in the numerator of @ 
define the positions of the renormalon singularities. The UV renormalons are 
given by the poles of F(— 1 + 2u) and are thus located at u = |, 0, — |, . . .0 They 
are due to the contributions of virtual momenta above fi in loop integrals. Note 
that there appear also IR renormalons (at u = 1,2, . . .), which are due to the 
contributions of virtual momenta below the soft scale k. Since, by construction, 
the effective theory and the full theory describe the same dynamics in the IR 
region, these renormalons reflect truly non-perturbative effects of QCD and are 
thus independent of the UV regularisation scheme. They are not related to the 
construction of the HQET and are therefore not the topic of our discussion. The 



*In this paper we define the running coupling constant in the so-called V scheme [^^, in 
which the one- loop counterterms differ from those of the MS scheme by an additive constant. 
Otherwise, the Borel transform has to be multiplied by a scheme-dependent factor e"*^". 

^The pole at u = is removed by renormalisation, however. 
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UV renormalon at m = 1/2 is the nearest singularity on the integration contour 
of the Borel integral in (H). An expansion of the Borel transform around this 
point leads to 

S;,\k,u) = -^-^^ + .... (6) 

The linear dependence on fi reflects the linear UV divergence of the inverse prop- 
agator in the effective theory, which in dimensional regularisation is not seen in 
perturbation theory. Using the above result together with (|^), we find that the 
renormalon ambiguity in the definition of the inverse propagator is of order Aqcd • 
Let us now compare these results to the one-loop calculation of the inverse 
propagator performed with a hard Pauli-Villars cut-off A. In this case one finds 
(for details see sect. |]) 

5i'W = -%^A + „.t{l + £^(i^lu^ + co„st.)}+0(a2). (7) 

We have given the result for an arbitrary covariant gauge. The Landau gauge 
corresponds to a = 0, whereas the Feynman gauge corresponds to a = 1. Note 
that the coefficient of the logarithmic term is the same in and (|^). The linearly 
UV divergent term in (|^) is in correspondence with the term proportional to fi 
in (P), which leads to a renormalon ambiguity of order Aqcd- Note that in the 
case of the cut-off regularisation the power divergence leads to an ambiguity of 
the same order, as explained above [see (|l|)]. 

This correspondence between UV renormalons and power divergences works 
in several cases. However, an apparent exception is provided by the quark matrix 
element of the kinetic energy operator {iD±)'^hj;, which appears at order l/mq 
in the expansion of the inverse propagator.^ This dimension-five operator has 
the same quantum numbers as the lower-dimensional operators (dimension- 
three) and iv-D (dimension-four). Hence, in a generic regularisation scheme 
with a dimensionful regulator, one expects a mixing with these operators, leading 
to quadratic and linear UV divergences. Indeed, these divergences appear in the 
lattice formulation of the HQET, where, at one-loop order, one findsQ [|16| 

S^^{k) = {k\K {iD^fK \k) = ^ + vk + .... (8) 

Here W and X are dimensionless constants, and the UV cut-off is provided by 
the inverse lattice spacing, A = 1/a. The puzzle that we wish to point out is that 

^Here hy denotes the velocity-dependent heavy quark field in the HQET, and D^,^ = 
— v-DVfj, contains the "spatial" components of the covariant derivative. For details see 
ref. and references therein. 

^ The calculation was originally performed for a heavy quark at rest. The generalisation to 
an arbitrary four-velocity is possible, although there are considerable subtleties in formulating 
the effective theory at non-zero velocity in Euclidean space 
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the quadratic divergence is not seen in some other cut-off regularisation schemes 
(such as Pauh-Villars) . Moreover, in the bubble approximation we do not find 
an UV renormalon pole corresponding to the quadratic divergence. The nearest 
UV singularity in the Borel transform of S^^l^{k) is located at -u = |, and an 
expansion around this point gives 

&,u) = ^-^-^v.k + .... (9) 

This behaviour corresponds to the linearly divergent term in (|^). However, we do 
not find a term of the form yU^/(u — 1), which would correspond to the quadratic 
divergence. We shall refer to this "missing" renormalon pole at m = 1 as the 
invisible renormalon. 

There are several possible explanations for this puzzle. The general argument 
given above is that there is a relation between the degree of divergence of the 
matrix elements of higher- dimensional operators and the position of the singulari- 
ties in their Borel transform in dimensional regularisation. The coefficients of the 
divergent terms and of the renormalon singularities, however, are not universal 
and may be zero in some regularisation schemes. Indeed, there is no reason why 
renormalon singularities or power divergences should be the same for all versions 
of the theory, since the ambiguities related to these renormalons (or powers di- 
vergences) are spurious and cancel in the predictions for physical quantities. It 
is therefore conceivable to find some regularisation schemes in which at one-loop 
order the coefficient of the quadratic divergence of the kinetic operator happens 
to be zero.0 However, it is surprising to us that we could find only a single regu- 
larisation scheme in which this coefficient does not vanish. Moreover, even if one 
accepts this point of view one still has to explain the absence of a renormalon 
pole at M = 1 in the Borel transform (P). A possibility would be that such a 
singularity appears only when one goes beyond the bubble approximation. 

On the other hand, one may wonder whether the appearance of a quadratic 
divergence in the lattice regularisation scheme is due to the breaking of some 
symmetry preserved in other regularisation schemes. In this case the invisible 
renormalon would be nothing more than a lattice artefact. Examples of such 
broken symmetries are Lorentz symmetry and, as a consequence, the so-called 
reparametrisation invariance of the HQET [^, ^ . We note that the breaking of 



reparametrisation invariance is responsible for the fact that the renormalisation 
constants of the operators h^iv ■ D hy and {iD±Yhy differ by a finite term of 
order as in the lattice version of the HQET [|16|, whereas they are the same in 
the continuum version of the theory |3l|. However, reparametrisation invariance 



does not forbid the mixing between the operators {iDj_yhy and h^, which 
is responsible for the quadratic divergence. 

If this is the quadratic divergence is Hkely to appear at two- loop order, since at 

one-loop order the operator (iD±)^h^ mixes with h^iv-Dh^, which itself mixes with h^. 



6 



This concludes our main discussion. In the following sections we present the 
details of the calculations, whose results were referred to above. 



3 Expansion of the inverse propagator 

Let us investigate the appearance of renormalons in the l/rriQ expansion of the 
inverse heavy quark propagator, which is the simplest Green function for a heavy 
quark. We shall generalise the analysis of Beneke and Braun by including 
terms of order l/mq. We work in dimensional regular isat ion and adopt the 
bubble approximation to investigate the singularities in the Borel plane. For 
technical details of the calculation the reader is referred to refs. |0, ^, |T2| . 
We start from the quark propagator in the full theory, 

S{p,m) = - ^— -, (10) 

— m — h[p, m) 

where m denotes the bare mass, and 

S(p, m) = m Si(p^, m^) + {jj — m) S2(p^, m^) (11) 

is the self-energy. We write the heavy quark momentum as p = rrigv + k, where 
rriQ is the expansion parameter of the HQET (in general mg ^ m), and consider 
Ei^2 = ^i,2{k,'m) as functions of the residual momentum k. Next we define a 
projected propagator Spik^mq) from the relation^ 

5p(/c,mQ) = — ^5(p,m) . (12) 

Our goal is to construct the heavy quark expansion of Sp^{k,mQ). Including 
terms of order I/itlq, we find 



Sp\k,mQ) = (l-S2)|(mo-m)+t;-fc + ^ (1-5/2) '} -Si m + 0(l/my , 

(13) 

where kj_ = k"^ — {v ■ k)"^, and 



'1 + Si - S2) m 
(1 - S2) niQ 



(14) 



The result simplifies if we restrict ourselves to the bubble approximation, which 
corresponds to the first term in an expansion in powers of l//3o. We can then use 
the fact that Si, S2, S and (mg — m) are of order I/Pq. At the same time we 

^The authors of ref. Q have projected the inverse propagator instead of the propagator 
itself. This procedure is not apphcable if one wants to include l/mg corrections in the HQET. 



7 



can substitute the bare mass m by the HQET expansion parameter mq in the 
expressions for Sj. This leads to 



T.i{k,mq) - Si(0,mQ) 
2mc 



mq 

+ 0(l//3o^ , (15) 



Sp^{k,mq) = -6m 

+ 1 - S2(/i;, mg)] {^v ■ k + 

where Sm = rripoie — mq, and 

mpoie = m {l + Si(0, mg)} + 0(1/ P^) (16) 

is the pole mass. The heavy quark expansion is consistent as long as the so-called 
residual mass 6m is a parameter of order Aqcd [p3 |. 

In order to see the appearance of renormalons, one has to consider the Borel 
transform of the inverse propagator. In the bubble approximation, explicit ex- 
pressions for Borel transforms of the pole mass and the functions Ei(/c, mq) have 
been derived by Beneke and Braun [Q. By expanding their results in powers of 
1/mq, we obtain 



Sp^{k, mq, u) = —6m{u) + 6{u) iv-k + 



k' 



2m, 



+ 



V ■ k 



+ 



Q 



[1 - u^) Ri + R2 



2m, 



Q 



[l-u^)Ri + (1 -2u) R2 



(17) 



+ 



{v ■ kf 
2mn 



Xiu)R,-3R2 \+Oil/Pll/ml) 



where 



6m{u 
mpo\e{u 

Ri 
R2 

X{u) 



rripoiciu) - mq6{u) 
6Cf 



(u) = m < 6{u) + 



f3o 



(1 - u) Ri 



2u 



+ RmM 



/mQ\ -2" r(n)r(l-2n) 



r(3 - u) 



-2v ■ fc\ -2" r(i - u) r(-i + 2u) 

Ji / T{2 + u) 

6 + 5u-2u^- 8u^ - Au^ 



2(1 + 2m) 



;i8) 



The scheme-dependent function Rm{u) is irrelevant to our discussion. It is con- 
venient to rewrite (|T7p in the form of a convolution of Borel transforms, which is 
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defined as 

u 

f{u) * g{u) = f ■ g{u) = J du f{u') g{u - u') . 



Then the result takes the form 

1 



(19) 



+ 



Ckin("^Q, U) * ^kin(^> W) + C(„.D)2 (mg, u) * S^^,^y {k, u) 

(20) 



2mQ 
+ 0{l/ml) 



in which the heavy quark expansion is explicitly realised, i.e. the dependence on 
the two scales ttlq and k is disentangled. The dependence on the heavy quark 
mass rriQ is contained in the functions 



ZQ\mQ,u) = S{u) + 



6Cf 



[l-u'')R,- — + Rz{u) 
2u 



C],m{mQ,u) = 6{u) , 
QCf 



C(^^.D)^{mQ,u) 



X(u) Ri 



3 

2u 



Rc{u) 



(21) 



which are the Borel transforms (in the bubble approximation) of the wave- 
function renormalisation factor Zq^ and the Wilson coefficients of the operators 
hv {iD±yhv (kinetic energy operator) and {iv-D)'^h^ in the effective Lagrangian 
of the HQET @, |3^. Note that = S{u) implies Ckin = 1, as required by 
reparametrisation invariance |31|. The dependence on the small scale k in (pO|) 
resides in the quantities 



Sj{k,u) 

'S'kin(^) 



V ■ k |5(m) + 
kl |5(u) + 
+ {v ■ kf 



6Cj 
6Cf 



i?2 + — - Rz{u) 

2u 



;i-2u) i?2 



2u 



Rz{u) 



3/^2 - — - Rc{u) 
2u 



(22) 



which are the Borel transforms (in the bubble approximation) of the quark matrix 
elements of the operators appearing in the effective Lagrangian of the HQET. 
Their calculation will be outlined in sect. Note that it is sufficient to keep 
terms of order unity in S(^.£))2 , since the corresponding Wilson coefficient 0(^0)'^ 
in (|2l|) starts at order I/Pq. The functions Rziu) and Rc{u) in (^) and (^) 
depend on the renormalisation scheme. They are irrelevant to our discussion. 
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It is instructive to trace how renormalons are introduced in the construction of 
the heavy quark expansion, and in which way they cancel between the coefficient 
functions and matrix elements. To start with, we note that even in the full 
theory the self-energy contains IR renormalons from the contributions of virtual 
momenta below the scale k. Their positions are determined by the factor r{l — u) 



contained in R2 in (18). Note that the leading IR renormalon at m = 1 is not 
forbidden by gauge invariance, since the self-energy is not gauge invariant and 
there exists a gauge- variant operator of dimension two [|2[. What is introduced 
in the process of constructing the heavy quark expansion are new renormalon 
singularities at half-integer values of u. In the coefficient functions there appear 
IR renormalons at positions u = |, 1, |, . . . determined by the factor r(l — 2m) 
contained in Ri. Likewise, in the matrix elements there appear UV renormalons 
at positions u = |,0, — determined by the factor r(— 1 + 2u) contained 
in i?2- Since these singularities are absent in the original theory, they must 
cancel, order by order in l/ruQ, in the effective theory. Consider now our results 
([l7|) and (|20|) to see how the IR and UV renormalon poles conspire. At order 
l/rriQ only the renormalons at u = 1/2 need to be considered. We observe that 
the IR renormalon in 5rn matches with the UV renormalon in S~g, and the IR 
renormalon in the coefficient Zq^ which multiplies matches with the UV 
renormalon in S^^l^. This becomes explicit if we expand (p!7|) around u = 1/2: 

bp {k,mQ,u) = -niQ 



8Cp ^ ' ^ ^2(l-2n) TTiQ 

3 yU 1 /i 



+ V ■ k 



A{l-2u)mQ 2{2u-l)v-k 



jv ■ kf 3 

2mQ 2{2u -\)vk "' ' ^ ^ 

where the ellipses represent terms that are regular at m = |. Notice the absence 
of an U V renormalon at m = | in the term proportional to k^ in S^^^ given in (p^ . 



Such a renormalon is forbidden, since it would give rise to a non-local behaviour 
of the form fikj_/v ■ k. On the other hand, nothing forbids an IR renormalon at 
M = 1 in (5m, which could conspire with an UV renormalon at m = 1 in the term 
proportional to (f ■ k)'^ in The fact that this does not appear is the puzzle 

mentioned in the introduction. 

Before we proceed, let us note that from an expansion of the expressions in 
( |2TD and (p2D around u = one can extract the one-loop expressions for the 

^°There is in fact a pole at u = 1 in S^^l^, but it has nothing to do with the UV region. 
Computing the matrix element of hy (iD±)^hy with a hard IR cut-off we find that this pole 
disappears. Hence, it is related to the region of momenta below k. This IR renormalon pole 
appears also in the full theory. 
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corresponding quantities. Keeping only logarithmic terms, we find 



Z^{m^) = l-- ^In — 



C(l.D)2(mQ) 



3(3 - g) Cpas mq 
2 n fi 



Sj{k) = v.k\l+^—^^ln^^ 



TT 



2 TT 

3(3 — a) Civets 



TT 



/i 

f ■ fc)^ In - 



-2v ■ k 



(24) 



where we have given the results for an arbitrary covariant gauge. The expression 
for the Wilson coefficient C(„.£))2 has been derived (in the Feynman gauge) in 



ref. |5|. 



4 Renormalons in HQET matrix elements 



In this section we show in more detail that the functions in (|2^) can be identified 
with the Borel transforms of quark matrix elements of operators in the effective 



Lagrangian of the HQET. At order I/ttlq, this Lagrangian reads |3i, 35 

1 



-^HQET = Kiv-Dh^ 



2mc 



9s J 



+ C^.,^K(r,,G^''K +0(l/m^), (25) 



with Ckin = 1 by reparametrisation invariance |3T|. The so-called chromo- 
magnetic operator a^^G^'^h^ plays no role for our discussion here, since its 
matrix element between quark states vanishes. We will now outline the calcula- 
tion of the matrix elements of the other operators between heavy quark states 
with velocity v and residual momentum fc, using the Landau gauge. 

The Borel transform of the inverse heavy quark propagator in the effective 
theory has been calculated by Beneke and Braun p[. The result is 



B T f 6C 

S^g{k) = {k\hviv-Dh^\k) v-k<.6{u)-\ — R2 



(26) 



After UV renormalisation, which amounts to removing the pole at m = con- 
tained in R2, this leads to the expression for Ses{k,u) given in (p^). 

Let us now turn to the calculation of the Borel transform of the matrix element 
of the kinetic energy operator {iD±)'^hy. The relevant diagrams are depicted 
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in fig. |T]. The shaded bubble represents the Borel transform of the resummed 
gluon propagator in the Landau gauge. In the bubble approximation, the Borel 
transform of any one-loop diagram is simply obtained by using this propagator 
instead of the usual one 0, [1^. We find that the contribution of the "vertex 
diagram" shown in fig. [I|(a) is 

6Cf 



{l-2u) ki-5{v-ky 



(27) 



The two "sail diagrams" depicted in fig. |l](b) and (c), which have a gluon attached 
to the operator insertion, each give a contribution 

6Cf 



R2 {v ■ ky 



(28) 



In dimensional regularisation the tadpole diagram of fig. |l|(d) vanishes. Adding 
the tree-level contribution to the above results, we obtain 



S^^{k) = {k\K{iD^yK\k) ^-5- kl6{u) + 



6Cf 



R2 



2u) k\ 



3 {v ■ ky] , 

(29) 

in agreement with the expression for S]^^{k^u) given in (^2]). 

Repeating the same calculation for the operator hy [iv ■ Dyh^, we find that 
the vertex diagram vanishes, and only the sail diagrams give a non-vanishing 
contribution. The result is 



S^^'ryrik) = {k\K{ivDyK\k) {v ky \6{u) + 



R, 



(30) 



which generalises the expression given in (p^ to order 1/(3q. 



5 Power divergences in HQET matrix elements 

In this section we present the results of several one-loop calculations of the HQET 
matrix elements S~^{k) and S^^{k) in regularisation schemes with a hard UV 
cut-off A. The aim is to study the correspondence between the UV renormalon 
poles encountered in the previous section and the power divergences associated 
with the use of a dimensionful regulator. For simplicity, we shall perform the 
calculations in gauges suitable for the regularisation method of choice. We note 
that the structure of the power divergences is gauge invariant. 



5.1 Pauli— Villars regularisation 

In the simplest version of the Pauli- Villars regularisation scheme, one substitutes 
for the gluon propagator in the Feynman gauge the expression 

= -W^u - -r^) = J^r^.o. ■ (31) 



2 q^-Xy g2(^2_p) 
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Computing the one-loop self-energy in the effective theory and adding the tree- 
level expression for the inverse propagator, we find (assuming v ■ k < 0) 

S- ' (*•) = 1, . 1 1 + ( 1,1 Z^Ji + arctari Vi^) } 

= _£|!lA + „.t{l + £^!l(,„z2^_l)}+0(l/A), (32) 

where x = \/{—v ■ k). The generalisation of this result to an arbitrary covariant 
gauge has been given in (0). The linear divergence corresponds to the UV renor- 
malon pole at m = | in (^). It is instructive to rewrite the result in a form that 
makes explicit the mixing of operators under renormalisation. We define 

S~^{k) = {k\hyiv-D hy\ k) 

= Z~Q {k\h^iv-Dhy\k)o + Z~^^^^ \ {k\h^h^\k)Q , (33) 

where {k\hyiv-D hy\k)Q = v ■ k and {k\hyhy\k)o = ^ are the tree-level quark 
matrix elements. From (|32|) we then read off the renormalisation constants 

(34) 

In order to regulate higher power divergences, it is necessary to introduce a 
more general form of the Pauli-Villars regularisation. Instead of (|3ll) , we shall 
use two subtractions and write 



^ . ^ . , 1 7 1,1 1 



(f 7 — 1 — 7 — 1 — 7A" 



Repeating the above calculation, we obtain 



(35) 



S-i{k) = v.k[l^^-I^ 



ln^ + ^/"\, (36) 



A 2(7 

H ^ {^\/ — 1 arctan Va;^ — 1 — y^x^ — 1 arctan y^a;^ — ij 



where = ^/yx = y/j X/{—v ■ k). This leads to 
^eff = 1 + ( In ^ + 



TT V A 2(7-1) 
V7 + 1' 



^1=-^^^- (37) 
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Eqs. ( p^D are recovered in the limit 7 cxd. 

The calculation of the matrix element of the kinetic energy operator is more 
complicated. Consider again the diagrams in fig. |I|, where now the shaded bubble 
represents the regulated propagator given in (|35D. We find that in the Feynman 
gauge only the vertex and the tadpole diagrams give non-vanishing contribu- 
tions. For dimensional reasons the tadpole diagram is quadratically divergent. 
Its contribution to S^^l^{k) is 

-^A^l^. (38) 
47r 7 - 1 

However, this quadratic divergence is exactly cancelled by the quadratic diver- 
gence of the vertex diagram. For the sum of all one-loop contributions we find 

'S'kin(^) — ^ ^ I ^ ^ (\/*^7 ~ ^ arctan y^x^ — 1 — — 1 arctan \/ x'^ — 

+ [kl - 3{v ■ kf 

-\ {\Jx'if — 1 arctan y^x^ — 1— 7^2^^ — 1 arctan Va^^ — 1 



, — 2t> ■ k In 7 

n ^ + — ^— (39) 

A 2(7-1) ^ ' 



1 / rr, 1 

Note that, in spite of the explicit factor of in the first term on the r.h.s., this 
term diverges only linearly with A. 

We are now ready to compute the mixing of the kinetic energy operator with 
lower-dimensional operators. We define 



S^l{k) = {k\K {iD^fK\k) 

+ ^d2^^.z) \'\k\h^iv-Dh^\k)Q^ ^d^^i {k\KK\k)o . (40) 



^kin {k\hv {iDi_fK I )o + Zj^l ^(^.D)2 {k\h^{iv-DfK\k)Q 



From an expansion of (^) in the limit of large A, we obtain 

1 CpCts /, —2v-k ln7 \ 

^.; = 1+^(1"—+^) 

7-1 _ SCpas -2v ■ k ^ ln7 

^Dl^iv-nr- ^ [^^ ^ 2(7-1) 



2 + 1 



^DLi = 0. (41) 



As mentioned above, there is no quadratic divergence in the sum of all diagrams, 
i.e. in Pauli-Villars regularisation there is no mixing of the kinetic energy op- 
erator with the operator h„. We have checked that this result holds true in 
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an arbitrary number of space-time dimensions. Note that the coefficient Zj^^^^ is 
related to the self-energy of the heavy quark by 



7-1 

^kin 



d{v ■ k) 



(42) 



This relation is a consequence of reparametrisation invariance |31|. It is satis- 



fied also for the expressions given in (H), which were obtained in dimensional 
regularisation. 



5.2 Momentum flow regularisation 

We shall now discuss another regularisation scheme in which a hard UV cut- 
off is used to compute the quark matrix elements in the HQET. It is based 
on the observation that the resummation of renormalon chains is equivalent to 
performing one- loop calculations with a running coupling constant [0]. Consider 
again the matrix element of the operator iv-D hy. Imagine performing the one- 
loop calculation of this matrix element with a running coupling constant. In the 
Landau gauge, the result can be written in the form 

00 , , — , 

S-iik)=vk[l + JdrwMr)^^^^], (43) 



where we assume that u = ~2v-k ^ Aqcd- It has been shown in ref. that the 
distribution function Wy-oiT) is related to the (unrenormalised) Borel transform 
of the inverse propagator. One finds 

\ - r ^ -u 3Ci. r(l -M)r(-1 + 2M) 





which can be inverted to give 

wMr) = ^ {(1 + 4t)'/' - 1 - 6r} . (45) 

The representation (^) offers a natural way to introduce a hard UV cut-off, 
since the product ^/TUJ has the interpretation of a physical scale. The distribution 
function Wv.d{t) controls the momentum flow in the one- loop self-energy diagram. 
For large values of r, this function falls off proportional to I/VT) so that the 
integral is linearly UV divergent. Introducing a hard UV cut-of^ ruv = A^/c<j^ 

^^In principle one should also introduce an IR cut-off, since the running coupling constant is 
not well-defined in the low-momentum region. This fact is irrelevant at one-loop order, however, 
where one can neglect the running of the coupling constant. 
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and expanding the result for A — > oo, we find for the renormahsation factors 
defined in ( |33D 

, 3Cf as /, -2t; ■ k 1 
+ — Vl'"^ 2 

(«) 



Note the similarity to the results (|34 ) obtained with a Pauli-Villars regulator. 

Let us now consider the quark matrix element of the kinetic energy operator. 
We define 



OO OO , — , 



TT 



(47) 

Using again the relation between the distribution functions and the Borel trans- 
forms given in (|2^), we obtain 

u)(vD)2{t) = -3idy.D{T). (48) 

For large values of r, the first function falls off proportional to 1/r, so that the 
corresponding integral is logarithmically UV divergent. After a straightforward 
calculation we find for the renormahsation factors defined in 

1 3Cf as ( -2v -k 1 



Z-l ,,._ = -l^:il(ln 



9Cf as ( , -1v ■ k 1 



'Dl-^{vDY 2 TT V A 2 

z 



^^^^i = 0. (49) 

Once again we find that the quadratic divergence, which is expected on dimen- 
sional grounds, is absent. In the momentum flow regularisation scheme, the 
absence of the quadratic divergence is directly linked to the absence of an UV 
renormalon pole at m = 1 in the Borel transform S^:^{k, u) in (|2^) , since the Borel 
transform determines the distribution function through an integral relation of the 
type (H. 

It is instructive to note the close similarity of the above results with (^). In 
both regularisation schemes and in an arbitrary covariant gauge, the results for 
the renormahsation factors can be written in the form 

„_i , , Cpas 1^3 - a -2v ■ k 



TT 



A 
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and 



where 



Z-J,^l = ~d, (50) 



^^Li = 0' (51) 



in Pauli-Villars regularisation (in the Feynman gauge), and 

cmf = — - 5 duF = 1 (53) 

in the momentum flow regularisation (in the Landau gauge). The coefficient c is 
gauge dependent, but the constant d is not. 

5.3 Lattice regularisation 

In this section we discuss the one-loop corrections to the propagator and to the 
kinetic energy operator in the lattice regularisation. All results reported in this 
section have been obtained in the static theory {v = 0) in Euclidean space. 

Among the possible lattice formulations of the HQET we choose the simplest 
one, in which the heavy quark propagator in a background gauge field U has the 
form 1^ 

Sl{x, y) = ^(x^ - y^) 5ix - y) V,{x\ y') , (54) 

where Vg{x^,y^) is the path-ordered exponential from (x, x"^) to {x,y^) along a 
path whose position in space is constant: 

P^(x^ y*) = Ul{x, x^ - a) Ul{x, x^-2a)... Ul{x, y^ + a) Ul{x, y^) . (55) 

At one-loop order, the calculation of the heavy quark self-energy in the Feynman 
gauge gives 



16, 37 



(5°(x, y)) = e{x' - y') 5{x -y){l + a^iXta'^ + Y)} + O(a^) , (56) 

where t = x"^ — y^, = QqI^t^ is the bare lattice coupling constant, and (. . .) 
denotes the average over the field configurations. The quantities X and Y are 
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given by 



X = - 



d^k 1 



Y 



4 J 27r2 A ' 

— 71 



with 



A = ^(1 — cos ki 



(57) 



(58) 



i=l 



The constant Y is logarithmically divergent. It contributes to the wave-function 
renormalisation and hence to the renormalisation of any local operator containing 
the heavy quark field. The quantity X is gauge invariant and has to be identified 
with the coefficient of the operator h^hy. Adopting the notation introduced in 
where now A = 1/a plays the role of the UV cut-off, one obtains 



-a.X 



(59) 



for the coefficient of the linearly divergent term. 

The corrections to the kinetic energy operator can be computed from 



Slix, y) = 9{x' - y') 5{x - y) M^'^ D\x, w') Vy\w\ y') . (60) 
For the operator we take 

(61) 



D\z, z') f{z) = -,Y. [U,{z) f{z + aj) + U]{z~ aj) f{z - aj) - 2/(. 
a i=i 

where j denotes the unit vector in the j-direction. At the tree level one finds 
{Sl{x,y))o = e{x'^-y^)6{x-y)— ^ ^ {5(x+aj-y)+5(x-aj-y)-25(f-y)} . 



W,4=y4 j = l 



At one-loop order one obtains [|1^] 

{Sl{x, y)) = 0(x^ - y^) S{x - y) as {Wta'^ + 2Xa-^] 



+ 



A + a, (Xta-' + Y)] {Sl{x,y))o. 



where X and Y are defined in (p7|) , A = 1 + | asW, and 

}d^k 1 



W = -Cf 



l + Ay - 1 



(62) 



(63) 



(64) 
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Using the notation introduced in (^0]), we find 

-asW . (65) 

Thus, in this case the mixing coefficient of {iD±)'^hy with hy does not van- 
ish. Moreover, since reparametrisation invariance is broken on the lattice, the 
renormahsation of the kinetic energy operator differs from the wave function 
renormahsation by a factor A 7^ 1. 



7-1 



6 Conclusions 

In the HQET the coefficients of higher-order terms in the expansion in inverse 
powers of the heavy quark mass are proportional to matrix elements of higher- 
dimensional operators. In general these matrix elements diverge in perturbation 
theory as powers of the UV cut-off, and using dimensional regularisation their 
perturbation series are not Borel summable due to the presence of UV renor- 
malon singularities in the Borel plane. The ambiguities due to the presence of 
these renormalons are cancelled by those due to IR renormalons in the Wilson 
coefficient functions of lower-dimensional operators. 

In this paper we have investigated in some detail the relation between power 
divergences and renormalons. This correspondence, which is based on dimen- 
sional arguments, appears to fail in one important case, that of the kinetic en- 
ergy operator. In a lattice regularisation h^j {iD_i_)'^hv mixes with hy with a 
quadratically divergent coefficient, as expected from the naive degree of diver- 
gence of the corresponding Feynman diagrams. However, at least to one-loop 
order in perturbation theory, this quadratic divergence is absent when using the 
Pauli-Villars and momentum ffow regularisations. Moreover, the correspond- 
ing UV renormalon is absent using dimensional regularisation. The absence of 
the corresponding IR renormalon in the pole mass had already been noted in 
ref . . These results suggest that there may be some symmetry, which prevents 
hv {iD±)'^hv from mixing with h^, and that this symmetry is broken with the 
lattice regularisation. An obvious candidate for such a symmetry is Lorentz (or 
Euclidean) invariance, which takes the form of reparametrisation invariance in 
the effective theory ||31], ^ . This symmetry does not seem to forbid such a mix- 



ing, however. Given the relevance of the kinetic energy operator to studies of the 



spectroscopy and inclusive decays of heavy hadrons [|T8|-[^, it is important to 
understand whether the absence of the power divergences and the corresponding 
renormalon is an accident of one-loop perturbation theory and the bubble ap- 
proximation, or whether it is a consequence of a more general principle. We hope 
that this paper will stimulate further investigation of this puzzle. 
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Figure 




(a) (b) (c) (d) 

Figure 1: Diagrams contributing to the quark matrix element of the 
kinetic energy operator (indicated by a square). 



23 



